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Abstract. We present results from infrared spectroscopic projects that aim to test the relation 
between the mass of a black hole, Mbh, and the velocity dispersion of the stars in its host-galaxy 
bulge. We demonstrate that near-infrared, high-resolution spectroscopy assisted by adaptive 
optics is key in populating the high-luminosity end of the relation. We show that the velocity 
dispersions of mid-infrared, high-ionization lines originating from gas in the narrow-line region 
of the active galactic nucleus follow the same relation. This result provides a way of inferring 
Mbh estimates for the cosmologically significant population of obscured, type-2 AGN that can 
be applicable to data from spectrographs on the next generation infrared telescopes. 
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1. Introduction 

The comparison of the star- formation history of galaxies with the accretion rate of black 
holes (BHs) at different redshifts is o f particular im portance for our understanding of their 
evolution (e.g jMarconi et al. 112003; iNetzer Il2009h . A challenge in this comparison is the 
determination of BH masses, Mbh, in intermediate and high-z galaxies, as well as in 
obscured active galactic nuclei (AGN). 

Fo r nearby sources, Mrh can be inferred from stellar orbits (e.g., iKormendv et al 



1998 : Genzel et al. l200fl Davies et al. 112006), gas or maser kin ematics (e.g., Mivoshi et al 



19951 ). and reverberation experiments (jPeterson et al. 1 120041) . Most of these techniques 



often fail for sources at z>0.5 due to angular resolution and sensitivity limitations of 
the existing instrumentation. At such redshifts, Mbh is often estimated indirectly, using 
the relation b etween the mass of a BH and the velocity dispers i on a of the stars in it s 



galaxy bulge (|Ferrarese fc Merritt I I2OO0I : iGebhardt et al. Il2000l : iGuftekin et"al~ll2009h 



The stellar a value is typically substituted by the velocity dispersion of the narrow-line- 
region (NLR) gas clouds, which is measured from the 5007A [O III] line assuming that 
the NLR kinematics are primarily determ ined by the gravitational potential of the bulge 



(IShields et al. Il2003t iGreene fc Hol f2005) 



Ho wever, the obscuration of the optical [O III] line can be non negligible (jKauffmann et al. 
20031 ). This can be a problem for NLR kinematic studies of obscure d, type-2 sources , 



which constitute a cosmologically significant AGN population (e.g., iGilli et al. 112007 ; 
Lacv et al. 1120071 ). We aimed to extend this technique in the MIR, where the NLR lines 
suffer little from obscuration. Moreover, the [Ne V] and [O IV] lines have a higher ion- 
ization potential than [O III] and are therefore likely to originate from clouds that are 
closer to the BH. 
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2. The existence of an Mbh — cr relation for the NLR gas emitting in 
the MIR. 

We analyzed Spitzer IRS and ISO SWS archival data of local AGN with Mbh mea- 
surements from reverberation experiments ( Peterson et al. Il2004l ). We detected resolved 
[Ne V] and [0 IV] line s at 14.32 /zm and 25 .89 /zm, respectively, in more than half of the 
sources in our sample (jPasvra et al. 120081) . The resolution-corrected a m easurements of 
both [Ne V] and [0 IV] follow impressively well the Mbh — cr relation ( Dasvra et al. I 
20081 : see Fi g.lHfor [Ne VI) with a scatter that i s comp arable to that of the stellar velocity 
dispersions ( Onken et al. 1 1200*41 iNelson et al. 1 12004 ) . supporting previous findings that 
the NLR gas is often gravitationally bound to the bulge. 

This result can have various applications for high-resolution MIR spectra from future 
IR observatories such as JWST and SPICA. Since the Mbh — cr relation holds for the 
NLR gas emitting in the MIR, it can be applied to distant and obscured sources with 
resolved MIR narrow lines to derive estimates of their BH masses. By comparing BH 
masses of type-1 and type-2 AGN at similar redshifts, it can also provide a testbed for 
the AGN unification model. 



The flux of [Ne V] or 



O IV] also correlates well with the optical 5100 A continuum flux 



( Schweitzer et al. 2006 ; Dasvra et al. 20081) for the observed AGN, indicating that it 



can be used as a pr oxy of the bolometric AGN luminosity as suming standard bolometric 



conversion factors (jElvis et al. Ill994t iMarconi et al. 1120041) . Therefore, high-resolution 



MIR spectroscopy can also provide a means to derive the Eddington accretion rates for 
distant and obscured AGN. 
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Figure 1. The relation between Mbh and the velocity dispersion o f a NLR emission line in the 
MIR, [Ne V], for reverberation- mapped AGN (jPasvra et al. |[2008f ). The solid line corresponds 
to the best-fit solution to the 
dashed line corresponds to the 



Ne V] velocity dispers ion values (excluding IC4329A) and the 
Trcmain e et al. I (|2002T ) relation. The stellar velocity dispersio n 
measurements of the same sources are plotted as st ars (lOnken et al. 112004 1; INelson et al. |[2004T ). 
Similar results are found for the [O IV] line width (jPasvra et al. Il2008f ). 



3. A way to improve the calibration of the local stellar Mbh — v 
relation. 

In order to determine the accretion rate of BHs at intermediate and high z using 
gas kinematics and the Mbh — cr relation, it is necessary to test whether the relation 
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remains identical with look-back time. Various studies in the literature suggest that an 
evolution of the r elation could be possible (e.g. JSalviander et al. 1120071 : iTreu et al. l[2007t 



Woo et al. 20081) . However, determining or quantifying such an evolution remains highly 



uncertain. The limitations are primarily originating from the selection biases of distant 
sources, and from the lack of a robust calibration of the AGN -Mbh — & relation in 
high-luminosity AGN in the local Universe. 

The sources that were initially used for the construction of the local AGN relation were 
primarily low-luminosity AGN. In QSOs, the absorption features that are traditionally 
used for the extraction of the stellar a, such as the Ca II triplet at 8498, 8542, and 8662A 
were undetected, heavily diluted by the bright AGN continuum. We demon strated that in 



such cases the stellar a can be measured from the .ff-band CO bandheads (jDasyra et al. 



20071 ) . The advantage of using NIR instea d of optical lines is that the AGN spectral 



energy distribution has a local minimum (jElvis et al. 1994) and the host galaxy has 
a local maximum close to 1.6 /um. The QSO s with Mbh and a measurements were 
located above the high-mass end of the relation ( Dasvra et al. | [2007: Wa tson et al. 
Fig. [2]). Differences in the kinematics of stellar populations alone could not account for 
the observed discrepancy. 
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Figure 2. The AGN Mbh sigma relation includi ng the host galaxies o f some 
reverberation-mapped QSOs. Adapted from lWatson et al. I (|2008T ) . 



It is possible that the calibration of the BH masses in QSOs requires modifications. 
Specifically, the statistically-determined value of the factor / that conver ts the virial mass 
enclosed in the broad line region to the BH mass (jOnkenet al. 1120041) could require a 



different prescription for AGN of high- and low- accretion rates (e.g. Collin et al. 20061 ) 



Alternatively, a steepening of the slope of the relation at the high lum inosity or mas s 



end is plausible, as it has been suggested based on [O III] kinematics ( Gaskell 20091) 



Given the uncertainties in the slope of the local relation, its redshift evolution should be 
examined with caution. 

To further study the behaviour of the high-luminosity end of the relation, we have 
acquired NIR integral field unit (IFU) datasets with SINFONI and NIFS (Grier et al. 
in preparati on) that are assisted by adaptive optics. An example, PG1426+015, was 
presented in IWatson et al. ( 20081 ). which shows a significant improvement in S/N ratio 
between the long-slit and the IFU spectra (Fig. [3]) . The success of the IFU data relies in 
spatially disentangling the AGN from the host galaxy spectra (Fig.|3|). 
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Figure 3. H-band spectra of PG 1426+015, obtained with the long-slit spectrograph ISAAC 
at the VLT (left panels; d ata fromlDasvra et alT1l2007f ) and the IFU NIFS on GEMINI north 
(right panels; data used in I Watson et al. II2008D . The stellar CO absorption features are more 
prominent in the IFU data despite the ~2 times shorter on-source integration time. 

4. Conclusions 

The combination of Spitzer IRS and ISO SWS spectroscopy enabled the study of 
NLR gas kinematics in the MIR for reverberation-mapped AGN. The stellar A/bh — c 
relation holds within the errors for NLR gas clouds that emit in the MIR as determined 
by the [Ne V] and [O IV] emission line widths. Therefore, resolved MIR narrow lines 
(in combination with the Mbh — cr relation) can be used to estimate the masses of BHs 
in obscured AGN. The measurement of stellar kinematics of QSO host galaxies can be 
achieved with NIR medium-to-high resolution long-slit and IFU spectroscopy. The hosts 
of local high-luminosity AGN have higher BH masses than those that the Mbh — o 
relation would predict for their velocity dispersions, similar to what has been found for 
high-z QSOs. This dictates the need for a better calibration of the relation prior to 
studying its possible evolution with z. 
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